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Abstract 

We examine the origin of neutrino masses and oscillations in the context of 
the six-dimensional standard model. The space-time symmetries of this model ex- 
plain proton stability and forbid Majorana neutrino masses. The consistency of the 
six-dimensional theory requires three right-handed neutrinos, and therefore Dirac 
neutrino masses are allowed. We employ the idea that the smallness of these masses 
is due to the propagation of the right-handed neutrinos in a seventh, warped di- 
mension. We argue that this class of theories is free of gravitational anomalies. 
Although an exponential hierarchy arises between the neutrino masses and the elec- 
troweak scale, we find that the mass hierarchy among the three neutrino masses is 
limited by higher-dimension operators. All current neutrino oscillation data, except 
for the LSND result, are naturally accommodated by our model. In the case of 
the solar neutrinos, the model leads to the large mixing angle, MSW solution. The 
mechanism employed, involving three right-handed neutrinos coupled to a scalar in 
an extra dimension, may explain the features of the neutrino spectrum in a more 
general class of theories that forbid Majorana masses. 



1 Standard Model in Six Dimensions 



The proposal that all the standard model fields access extra spatial dimensions above 
some energy scale ( "universal extra dimensions" ) has received considerable attention dur- 
ing the past year. Precision electroweak measurements require only that the compacti- 
fication scale of universal extra dimensions be above a few hundred GeV, opening up a 
potentially rich set of signatures, both in additional precision measurements p|, Q and in 
collider searches H, 

An especially attractive possibility is that there exist two universal extra dimensions. 
The six-dimensional standard model is chiral, and the constraints from Lorentz invariance 
and anomaly cancellation have remarkable consequences. The quarks {Q,U,V) and lep- 
tons {C,S) are four- component Weyl fermions of definite chirality, labeled by + and — . 
The cancellation of irreducible gauge anomalies imposes one of the following two chiral- 
ity assignments consistent with Lorentz invariant Yukawa couplings: Q+,U^,V^, C^,S±, 
where generational indices are implicit |Q. The reducible gauge anomalies can be can- 
celed via the Green-Schwarz mechanism as discussed in |^ |[ . Gravitational anomaly 
cancellation requires that each generation include a gauge singlet fermion A/± with six- 
dimensional chirality opposite to that of the lepton doublet 0. In addition, the six- 
dimensional standard model is the only known theory that constrains the number of 
fermion generations to be = 3 mod 3, based on the global anomaly cancellation condi- 
tion [§. 

The two universal extra dimensions have to be compactified on an orbifold, so that 
each of the six-dimensional chiral fermions gives in the effective four-dimensional theory 
either a left- or a right-handed zero-mode fermion. The simplest orbifold compactifications 
are either the square T"^ /Z2 or T"^ /Z^ orbifolds. 

An intriguing feature of the six-dimensional standard model is that the combination of 
its Lorentz and gauge symmetries can lead to a sufficient conservation of baryon number, 
even with the scale of baryon-number violating physics as low as the TeV range |Q . For 
the T"^ /Z^ orbifold, a Z^ subgroup of the six- dimensional Lorentz symmetry is exactly 
preserved. In the case of the square T"^ /Z2 orbifold, the same is true provided the two 
orbifold fixed points that are exchanged by a 90° rotation in the compactified (transverse) 
dimensions are physically indistinguishable. The Zg symmetry requires that the baryon 
and lepton numbers, AB and AL, of any operator in the low-energy four-dimensional 
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Lagrangian obey the selection ruleQ 

3AE + AL = mod 8. (1.1) 

As a result, the proton is very long lived (all AB = 1 transitions are governed by very 
high-dimension operators, and are therefore strongly suppressed), while neutron-anti- 
neutron oscillations {AB = 2, AL = 0) are forbidden. In the lepton sector, there are no 
neutrino Majorana masses^, and more generally neutrino-less double beta decays {AB = 
0, AL = 2) are forbidden. The absence of Majorana masses follows from the properties 
of the gamma matrices in six dimensions, namely that the charge conjugation operator 
does not flip the chirality. 

In this paper we study the implications for neutrino physics of the six-dimensional 
standard model. The mass matrix for the zero-mode neutrinos is induced dominantly by 
the following dimension-seven Yukawa terms in the six- dimensional Lagrangian: 

-T_\%Mlta2n* + h.c., (1.2) 

where label the generations, Ti is the six-dimensional Higgs doublet, and where we 
have taken the six-dimensional chirality of £ to be —. The ensuing Dirac masses of 
the three neutrino flavors can accommodate the neutrino oscillation data and all other 
experimental constraints, with the exception of the LSND result |TD[. It is nevertheless 
difficult to explain why the eigenvalues of the Yukawa matrix Xj^ are extremely small. 

Since the standard model in extra dimensions is an effective theory, breaking down at 
some scale Mg in the TeV range, it is natural to expect gravity to be strongly coupled there 
as well. A structure that accommodates the observed weakness of the gravitational inter- 
action should then be added to the universal extra dimensions. The simplest possibilities 



are that either some number of additional fiat dimensions [11 1 or one additional warped 



dimension IIT^II are transverse to the universal ones and are not accessible to standard 



model fields. Each of these alternatives also provides a possible mechanism for explaining 
small but finite Dirac neutrino masses, as first proposed in [l^ and |l^, respectively, by 
letting the gauge singlet fermions propagate in these extra dimension(s). 

^The cancellation of anomalies via the Green-Schwarz mechanism requires a (four-dimensional) scalar 
field that transforms nontrivially under the Zs- Operators that involve this field can be induced by four- 



dimensional instanton effects and could result in a violation of the selection rule Eq. (1.1). We expect 
these effects to be negligible. We thank E. Poppitz for discussions on this point. 

^The fact that the Zs symmetry forbids Majorana masses was not taken into account in Ref. 0|. We 
note that even if the two universal extra dimensions were compactified on an arbitrary T"^ /Z2 orbifold, 
an exact Z4 symmetry would still have prevented any Majorana mass. 
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We concentrate here on the possibihty that the singlet fermions, along with gravity, 
propagate in a single additional (seventh) warped dimension, with the standard model 
fields confined to a 5-brane. We adapt and generalize the five-dimensional model of 



Ref. and examine its consequences for neutrino masses and mixing angles. In Sec- 
tion 2 we present a simple effective theory, involving a single scalar field in the seven- 
dimensional bulk, that couples to the M fields. We then discuss, in Section 3, the global 
gravitational anomalies in this context and argue that there is no additional constraint 
on the number of M fields (contrary to claims made in the literature |T^, |15|). 

We examine the neutrino zero-modes in Section 4, where we derive their mass matrix 
in terms of the wave functions of the M zero-modes at the standard model 5-brane. In 
Section 5 we derive the shape of the scalar vacuum expectation value (VEV), and then 
find the profiles of the M zero-modes. The effect of the scalar VEV is to concentrate the 
M zero-modes near the brane opposite to the standard model one, so that the resultant 
four-dimensional neutrino mass matrix is exponentially suppressed relative to the weak 
scale. We then show that the mass hierarchy between the different flavors is limited by 
the presence of higher-dimension operators, which have a flavor mixing effect. This is an 
important result, especially in view of the often stated existence of a large exponential 
hierarchy induced between the masses of different flavors, whenever the two chiralities are 
localized at separate branes 0, [16 |. 

The implication, discussed in Section 6, is that the ratio between the mass scales as- 
sociated with the atmospheric and solar neutrinos is not expected to be larger than one 
to two orders of magnitude. Therefore, the large mixing angle MSW solution to the solar 
neutrino problem is a consequence of this model. We discuss the neutrino mixing angles, 
as well as mass eigenvalues, and show that values compatible with all current neutrino 
oscillation data, except for the LSND result, emerge naturally. The energy scales associ- 
ated with the seventh dimension, as well as the universal six dimensions, are such that 
the singlet-neutrino Kaluza-Klein (KK) modes are too heavy to play a direct role in the 
observed neutrino oscillations. Similarly, the constraints from astrophysics [|17| or cosmol- 
ogy [0 on the mass of the KK neutrinos are not relevant here. In Section 7, we discuss 
the relations among various parameters in the model and draw some conclusions about 
the expected mass scales that characterize it. In Section 8, we summarize the essential 
ingredients of our model leading to a viable neutrino mass spectrum, and emphasize that 
the mechanism applies to a more general class of higher-dimensional models. 
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2 A Warped Seventh Dimension 



The seven-dimensional gauge-singlet fermions are Dirac spinors with eight components, 
denoted by A/"*(x*^), where i = 1,2,3 labels the generations. The spacetime coordinates, 
X with M = 0, 1, ...,6, are labeled as follows: for the ordinary spacetime, 

x^, for the two additional universal dimensions, and x^ = z for the dimension inaccessi- 
ble to the standard model fields. We use the following conventions: capitals M, N . . . (from 
the middle of the alphabet) denote the seven coordinate indices in a curved background, 
while capitals A,B, . . . (from the beginning of the alphabet) denote the corresponding lo- 
cal Lorentz indices. We also use lower case greek letters a, n, ... to refer to the coordinate 
indices, and lower case latin letters, a,m, . . ., to refer to the Lorentz indices along the flat 
universal dimensions. 

The [T'^/Z^j X [S^/Z2) orbifold compactification projects out the unwanted zero- 
modes, and restricts the coordinates to < x^,x^ < ttRu and < z < Tir^.. The 
six-dimensional standard model fields are localized at z = vrrc, while the gauge-singlet 
fields propagate in the whole bulk. 

The most general metric consistent with six-dimensional Poincare invariance is diag- 
onal, and warped in the z direction. However, the compactification of the two universal 
extra dimensions on the T'^ /Zi^ orbifold breaks six- dimensional Poincare invariance, and 
in general leads to a warp factor for different from the warp factor for the familiar 

uncompactified dimensions. For example, we expect contributions to the stress-energy 
tensor, due to the Casimir energy of bulk fields, that do not respect the six-dimensional 
Lorentz invariance. However, we will assume for simplicity that these differences can be 
neglected.0 Our main conclusions do not change if we allow for the more general possibil- 
ity that the warp factors for the uncompactified and compactified universal dimensions 
are different. Therefore, we consider a diagonal metric Gmn that is warped in the z 
direction, corresponding to a line element 



where /x, = 0,1,. ..,5, and rj = diag(+l, — 1, — 1) is the six-dimensional Minkowski 

complete solution that incorporates the gravitational backreaction of Casimir energies or other 
effects arising from the compactification of the two universal dimensions would involve the specification 
of a radius stabilization mechanism. We leave such a study for future work. 



Gmn dx dx 



.N 
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metric. 

Starting in Section 5 we will take the warp factor to have the form 

w{z) = e^^^'--''^ , (2.2) 

which is a good approximation whenever the dominant contribution to the bulk stress- 
energy tensor is due to a bulk cosmological constant. This normalization is chosen so 
that w{T{r^ = 1, which facilitates the physical interpretation at the standard model 
brane, located at z = tcVc. In particular, this choice implies that the coordinate radius 
Ru is the proper radius of the universal extra dimensions as measured by standard model 
probes. For the AdS metric defined by Eqs. ( |2.1| ) and ( |2.2| ), the Riemann curvature tensor 
is Rxpau = -k'^iOapOxu " OupOxa), usiug the sign conventions of |19 . 



So far we have introduced three mass parameters: the inverse coordinate radius 1/Ru 
of the universal extra dimensions (associated with T^/Zi), the inverse radius l/vc of the 
dimension accessible only to neutral fields, and k. They are all taken to be below the 
fundamental seven-dimensional mass scale M*, which, as will be discussed in Section 7, is 
in the TeV range. We will see in Section 4 that, with the normalization w(7rrc) = 1, the 
mass scale for the standard model KK modes is set by R~^, and thus satisfies the bound 
Ru^ ~ 0.5 TeV, imposed by the electroweak data The radius of the dimension 
accessible only to neutral fields is rather loosely constrained by searches for new long-range 
forces. 

In addition to the A/'*(x*^) fermions and the graviton, other fields that are singlets 
under the standard model gauge group could be present in the warped extra dimension. 
To describe naturally small neutrino masses it is sufficient to include a single real scalar 
ip with the dynamics described in the framework of effective field theory. This scalar 
is thus an effective degree of freedom, and could well represent a composite structure, 
with the compositeness becoming evident at scales of order M^, and above. The seven- 
dimensional (effective) action, invariant under both general coordinate and local Lorentz 
transformations, is then given by 



/ 



d'^x (^^'T^e/''DM^f' - h.c.) + ^G^'^'dM^dj,^ - V^^M 

+ 5{z - Ts:rc) Csu\ , (2.3) 



where the first two terms are kinetic terms in the warped spacetime and the last two terms 
describe the bulk interactions of the A/" and (f fields, and the six-dimensional standard 



5 



model. Here e/^ is the inverse vielbein, G is the determinant of the seven-dimensional 
metric, with \/G = w^{z), and g is the determinant of the six- dimensional induced metric, 
with ^/—g = w^^nrc). The T"^ are the anti-commuting matrices in seven-dimensional 
Minkowski space: the gamma matrices of six-dimensional Minkowski space along with 
= iTj, where = T^...T^ defines six- dimensional chirality via A/± = ^(1 ± Ty)^/. The 
fermion covariant derivative in Eq. (p.3|), associated with the diagonal metric Gmn, is 



r w-\z)[da + i6\TaT7{dw/dz)/2]J\f , A = a = 0, 1, 5 
e/'DM^/ = (2.4) 
[ dAf/dz , A = 6 . 

The bulk interactions preserve the orbifold Z2 symmetry, defined such that A/1 and 
are odd, while Af^ are even. They may be organized into a tower of operators of increasing 
mass dimension: 

y.,-^ = - A - iKf'- + + (-^f - ifk^"^-^) + (2,5) 

where A is a bulk cosmological constant that needs to be fine-tuned in order to keep flat 
the four dimensional sections. At the classical level this involves tuning A against possible 
brane tension terms as well as the vacuum energy stored in the ip VEV. The parameter 
\^p is real, and h are hermitian matrices. All are dimensionless. The mass-square in the 
second term of V^^,Ar is chosen to satisfy > 0, so that ip has a nonzero VEV. Both M<^ 
and k are taken to be well below to justify the use of effective field theory for exploring 
the vacuum properties of (p. By a flavor transformation it is possible to diagonalize the 
first term inside the parenthesis, so that we can use a basis where 

hij = hiSij , (2.6) 

with hi real and positive. The other terms involving A/"* are in general flavor non-diagonal. 

The six-dimensional standard model Lagrangian, £sm, localized at z = irrc, includes 
the kinetic terms for the lepton and Higgs doublets and the Yukawa interactions of A^*: 

^SM 3 i£^_T''ea''D^Cl + g^^^D^H^D^H - i^^Z^Ulia^l-C + h.c. ) , (2.7) 

where the induced (inverse) metric and vielbein at the standard model 5-brane are given 

by 
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e " 



w-\T:r,)5a" , a = 0,1,..., 5. (2.8) 
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In Eq. ( p.7|) , Da are the gauge covariant derivatives, and the Yukawa couphngs are again 
dimensionless. Note that the four-component field £_ has mass dimension +5/2 while the 
A/+ field (also four-component), being defined in seven dimensions, has mass dimension 
+3. 

Before proceeding with the analysis of the neutrino masses, we discuss the consistency 
of the seven-dimensional theory. 

3 Gravitational anomalies 

We next show that the seven-dimensional model described in the previous section is 
anomaly free. The reader interested mostly in neutrino phenomenology may wish to 
move directly to Section 4. 

The seventh dimension is compactified on a IZ2 orbifold and the six-dimensional 
standard model is localized on a 5-brane at one of the two fixed points, while the three 
singlet neutrino fields propagate in the bulk. It was shown in [^] that if all fields were six- 
dimensional, the resulting theory would be free of gauge and gravitational anomalies, both 
local and global. Letting the neutrino fields propagate in a seventh dimension amounts 
to adding three infinite towers of KK fields to this theory. Since all gauge fields are 
localized at the orbifold fixed points, allowing the singlets to propagate in more dimensions 
cannot introduce any gauge anomalies. Gravity, however, propagates in the bulk and one 
must consider whether all gravitational anomalies cancel. When coupling fermions to 
gravity, there can be two types of anomalies: those associated with general coordinate 
transformations and those associated with local Lorentz transformations. For each of 
these cases one must distinguish between local and global anomalies.^ We analyze first 
the case of local gravitational anomalies. After showing that there are none we turn to 
the more subtle issue of global gravitational anomalies. 



^Note that the word "local" can have two different meanings. In the context of "local Lorentz" 
transformations it means that the Lorentz group has been gauged, the standard usage in gauge theories. 
There can be, however, local Lorentz transformations that are continuously connected to the identity as 
well as local Lorentz transformations that are not. It is customary to refer to the transformations of the 
first kind as local and to those of the second kind as global. We use italic fonts whenever we want to 
emphasize the distinction between the transformations that are continuously connected to the identity 
and those that are not. 
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3.1 Local gravitational anomalies 



A noninvariance of the effective fermion action under local Lorentz transformations would 
imply that the corresponding stress-energy tensor is not symmetric. This would be 
incompatible with general covariance and the conservation law VmT^'^ = 0. Thus, in 
the presence of local Lorentz anomalies, either general covariance is broken or T'^'^'^ is not 
conserved. Anomalies associated with general coordinate transformations, on the other 
hand, lead directly to VmT^'^^ ^ 0. In either case, the theory that results when gravity 
becomes dynamical is inconsistent, and it is necessary to ensure that all gravitational 
anomalies cancel. However, the conditions derived from the requirement of anomaly 
cancellation for both kinds of transformations are not independent. At least in the case 
of local anomalies, it is possible to shift the anomalies of one kind into anomalies of the 
other kind by adding suitable local terms to the vacuum functional [20, 21 1. Thus, we 
may concentrate only on general coordinate transformations. 

In the case of gauge theories in three or five dimensions, any local gauge noninvariance, 
which is necessarily localized at orbifold boundary points, can always be cancelled by a 



bulk Chern-Simons term ||2^, ^ , provided the anomalies in the lower dimensional effective 
theory vanish.^ We now argue that this is also the case for local general coordinate 
anomalies in seven dimensions.^ We follow the argument given in for the spin-1 
case. The idea is to calculate the one-loop contributions to the covariant divergence of 
the seven- dimensional stress-energy tensor in the six- dimensional effective theory. If we 
regularize in such a way as to produce the covariant form of the anomaly, it is possible to 
perform the calculation in any convenient gauge. 

For the analysis of anomalies, it is sufficient to consider small fluctuations about a flat 
background 

ds"^ = [rjafj + hai3{x, z)] dx°'dx^ — dz^ , (3.1) 

where h^^ ^ 1. In Eq. ( p.lD we took advantage of the gauge freedomQto set G^z = and 
Gzz = —1- We also choose the vielbein as follows: Caa = ijlaa + \haa), ^az = ^za = and 

^The mathematical relation between Chern-Simons forms in odd dimensions and anomalies in even 
dimensions was discussed in |2^, 25 . 



^In five dimensions there are no local gravitational anomalies: the triangle diagrams always vanish. 

^The invariance of the line element Eq. (^]^) under the reflection z —z requires Gap and Gzz to 
be even, while G^z should be odd. For consistency, the infinitesimal parameters of a general coordinate 
transformation (C^) should be even (odd). Although these boundary conditions do not allow the zero 
mode of Gzz to be gauged away, this "radion" mode has vectorlike couplings in our theory so that it does 
not contribute to the anomaly, and we do not include it here. 
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= 1- The fact that we take a symmetric vielbein means that the stress-energy tensor 
is symmetric. (In this gauge the absence of local Lorentz anomalies is explicit.) 

The action for a fermion in the background Eq. (|3.1| ) becomes 



S 



d^xe^r%^DM^ + h.c. 



where 



T 



MN 



(3.2) 



(3.3) 



and all components with an index along the seventh dimension vanish. Here we use the 
notation '^d^ = '^d'^ — ((9\E')\E', and in the second line of Eq. (|3l2|) as well as in Eq. ( ^I3D 
it is understood that all indices are raised and lowered with the Minkowski metric rja/B- 

We now expand the fermion fields in KK modes: 



(3.4) 



where 4'±{^j ^) = ^ ^7)'4'{x, 2;), with = —iT^. The KK wavefunctions ^^(-z), which 



can be taken real, are solutions to 



(3.5) 



where ^^(-2) and ^„ (z) satisfy Neumann and Dirichlet boundary conditions, respectively, 
and are normalized as 

pTrrc 

dz^^iz)^^,iz)=6nn' . (3.6) 



(3.7) 



The result of replacing the mode expansion Eq. ( p.4| ) in the action Eq. (|3.2| ) is 



where 

rp{n,n') 1 

and 



(«') 



dzii{z)ii,{z)h^P{x,z) . (3.8) 

The action Eq. (|3.7|) corresponds to the six-dimensional theory of an infinite number of 
fermion fields that couple (chirally) to background fields h'^^ with standard gravitational 
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couplings. Note that when ha/s is 2;- independent, the resulting couplings are vectorlike 
(except for those of the zero-mode fermion) due to the normalization condition Eq. ( |3.6| ), 
which is the same for both chiralities. 

One can calculate now the relevant square diagram with one insertion of the 
operator Eq. (|3.3|) , with the seven- dimensional fermions replaced by their KK mode ex- 
pansions. Performing then the same manipulations as in p3|, and adding the contribution 
of a brane fermion (in our model these are the electrically neutral component of the £L, 
while the corresponding bulk fermions, labeled generically in this section by are the 
A/"*) to compensate for the zero-mode projected out by the orbifold boundary conditions, 
one can finally write 

/ ci^CMV^T^^^ = i [6{z - vrrj - 6{z)] ^ / {^Tr [v,R] Tr [r'] + ^Tr [v,R'] } . 

(3.9) 

Here we used a compact differential form notation: R'^jj = \R°'(3^y dx^ A dx'^ , where 
R°'i3^i,{x, z) is the Riemann tensor calculated from the background Eq. ( |3.1| ), but with 
the indices running only from to 5. The traces are taken over the indices that are 
not saturated by differentials. Also, (f^)"^ = dpC^, where the C^{x,z) can be thought of 
as the infinitesimal parameters of a general coordinate transformation. Due to general 
covariance, the result Eq. holds in any gauge. In addition, we are allowed to replace 
the Riemann tensor in Eq. ( p.9|) by its exact, nonlinear expression, so that the final result 
holds in an arbitrary background. 

The covariant anomaly given in Eq. ( |3.9| ) does not satisfy the Wess-Zumino consis- 
tency conditions and therefore cannot be obtained from the general coordinate variation 
of a functional of the metric. There is a standard procedure to obtain the consistent 
anomaly by adding local terms to the stress-energy tensor |2^. The resulting anomaly 
Ql{v^, r, R)[6{z — nrc) — 6{z)]/2, where F^/j = V^p^dx^ is the connection 1-form, is related 
to the variation of a 7-form Q-j{T,R) that can be added to the seven-dimensional action 
by 

5,^ j Q7(r, R)= j dQliv^, r, R) . (3.10) 

The Chern-Simons secondary characteristic class Qi{r,R) involves traces over all seven 
dimensional indices. | However, only the six-dimensional components contribute to the 
right hand side of Eq. (|3.10|) due to the orbifold boundary conditions on the metric and 

^Q7iT,R) satisfies dQriV^R) = -(2^ {^r [R"^])^ + ^Tr [R^]^ where R is the seven- 
dimensional Riemann curvature two- form. For a pedagogical exposition see . 
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on the infinitesimal parameters C,^ . 

The modified (consistent) form of Eq. ( p. 91 ) matches precisely with Eq. ( tj.l(J| ). We note 



that Qj{V,R) is odd under parity (defined as refiection through the z = hyperplane). 
Therefore, if we define the orbifold theory by starting from a compactification on the 
circle S^, the coefficient of the Chern-Simons term must change sign when crossing z = 
(and z = nrc), so that the theory is invariant under the refiection that is used in the 
orbifold projection. Due to this discontinuity, the gauge variation of such a term gives 
rise to delta-function singularities as in Eq. (|3.9| ). Alternatively, we can think of the 
/Z2 orbifold as a compactification on an interval (the half circle) with certain boundary 
conditions imposed at the end points. In this picture the coefficient of the Chern-Simons 
term is constant and the compensating anomaly comes from the boundary contributions. 
In either picture, it is clearly possible to cancel the noninvariance of the original fermion 
effective action by including the seven- dimensional Chern-Simons form Qj{T,R). We 
therefore assume that this Chern-Simons term is present so that the vacuum functional 
is invariant under local coordinate transformations. 

3.2 Global gravitational anomalies 

There still remains the question of global general coordinate/local Lorentz transforma- 
tions If there are diffeomorphisms not continuously connected to the identity, 
the previous analysis is not enough to ensure that the theory is invariant under such 
transformations. We phrase the following analysis in terms of general coordinate trans- 
formations, but the same arguments apply for the case of local Lorentz transformations. 

If W{G) denotes the fermion determinant in the presence of a background metric G, 
we have in general 

w(np\ 

e'^w , (3.11) 



W{G) 

where C denotes the metric obtained from G under a representative p of one of the 
disconnected classes of diffeomorphisms, and 5[p] is a phase that depends on the class to 
which p belongs.^ We are specifically interested in the diffeomorphism classes of = 

^We note that the group of disconnected diffeomorphisms of the n-sphere S*" is trivial for n < 5. This 
follows from the absence of "exotic" (n+ l)-spheres for n + 1 < 6 |^. There are therefore no constraints 
on the number of bulk neutrinos in the popular five-dimensional jZi orbifolds from global anomaly 
considerations. For higher dimensional theories the situation is not as straightforward, as indicated by 
the fact that there are 28 disconnected components on and two on S"^ . See for some other higher 
dimensional cases. 
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S'^ X [T"^ /Z4) X [S^ /Z2). The key observation is that S'^ /Z2 is diffeomorphic to an interval 
and any diffeomorphism of M. onto M. can be continuously deformed into one which is 
trivial in the interval jZ^- 



x'^ = p(x) z' = z, (3.12) 

where a;" denote the coordinates in 5*^ x {T'^/Z^). We can therefore restrict attention 
to coordinate transformations of the type (|3.12| ). It is then convenient to perform a KK 
decomposition for all fields, including the background metric, and analyze the resulting 
six- dimensional theory. Regarding the background, we note that under the limited class 
of diffeomorphisms ( |3.12|) the affine connection transforms as 

_ dx'^ dx^ dx'' dx'^ O'^xP 

Q^p Q^la Q^,p + Q^p dx"^dx'P ' ^ ' ' 

while all other components transform simply as tensors (the second term in Eq. ( p.l3[ ) 
would vanish if any of the indices a, /?, A were along the seventh dimension parameterized 
by z). Furthermore, if we perform a suitable KK mode expansion, only the zero mode of 
is affected by the second term in Eq. ( p.l3|) ; all other KK modes are true tensors under 
(|3.12|) . (When referring to the gravitational background field, a zero mode is defined to 
be independent of the higher- dimensional coordinate.) 

We next argue that if there are no global anomalies in a purely zero-mode gravitational 
background, then there are no global anomalies even in the presence of the higher gravity 
KK modes. The reason for this is that the group of disconnected diffeomorphisms is 
finite, at least for the case of the n-sphere 5*" |3^. It follows that for any element p in 
the group there exists a (smallest) integer such that is the identity element, and 
therefore the phase in Eq. ( p.ll|) associated with p must be an integer multiple of 2ti/N . 
If this phase vanishes when the higher gravity KK modes are turned off, and we turn 
them on smoothly, the phase must remain zero, unless it changes discontinuously. Since 
the higher gravity KK modes are just like background "matter" fields in the appropriate 
tensor representation of (|3.12| ), we consider this very unlikely. 

Now, in a zero-mode gravity background, the theory in question is just the six- 
dimensional standard model with the addition of three infinite towers of massive neutrino 
KK modes, which have vectorlike couplings to the background gravity field. We also note 
that the Chern-Simons term that is needed to cancel the local anomalies is invariant under 
( |3.12| ) when the higher gravity KK modes are turned off. Therefore, the fermion effective 
action in a zero-mode gravity background is invariant under general coordinate transfor- 
mations. From the argument given in the previous paragraph it follows that there are no 
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global anomalies in an arbitrary gravitational background. It is worth pointing out that 
in the presence of the higher gravity KK modes, the Chern-Simons term QjiT, R) is not 
invariant even under the restricted class of diffeomorphisms (|3.12| ). This noninvariance 
must be canceled by the rest of the terms involving the higher gravity KK modes. 

We conclude that adding a seventh dimension compactified on /Z2 to the six-dimen- 
sional standard model and letting the neutrinos propagate in the bulk, introduces neither 
local nor global gravitational anomalies. Turning the argument around, we can say that 
the consistency constraints on the number of neutrinos in the seven- dimensional model 
are the same as in the six- dimensional standard model analyzed in 0, namely it is nec- 
essary to include one singlet neutrino per generation. From the point of view of anomaly 
cancellation it is immaterial whether these neutrinos are bulk or brane fields. 



4 Neutrino Masses 

We now return to the action of Eq. (|2.3| ), which leads to Dirac neutrino masses. In order 
to study the implications for neutrino oscillations it is sufficient to analyze the zero-mode 
spectrum. 

The zero modes of A/"! are projected out by the orbifold boundary conditions. The 
KK decomposition along the warped dimension that includes the zero-mode M fields is 
given by 

A^Kx", z) = ^ AC^")(x°)C"(^) , (4.1) 



V^^c n=0 

where the index a = 0, 1, 5 labels the universal dimensions. The Q{z) form a complete 
set of orthogonal (dimensionless) functions on the [0, vrrc] interval, satisfying Neumann 
boundary conditions appropriate for even fields. They are chosen to obey the ortho- 
normality conditions 

dzw\z)ir{z)S:nz)=5nn' , (4.2) 



1 

TTVc Jo 

which ensure the canonical normalization of the six-dimensional kinetic terms for A/'+"^ 

We now adopt a warp factor chosen to be unity at the standard model brane, as 
in Eq. ( p.2|) . All kinetic terms for the standard model fields are then automatically 
canonically normalized. Keeping the zero-modes with respect to the seventh dimension 



only, and integrating over z, the six-dimensional effective Lagrangian [see Eq. (2^)] is 
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^M^^r ATf ^^a^H* + h.c. 1 + ... , (4.3) 

where a = 0,1,..., 5. Note that we do not need to distinguish between coordinate and 
Lorentz indices anymore, and all the indices are raised and lowered with the flat metric 



ria/B- Eq. ( [4.3|) shows that the mass scale for the standard model KK modes is set by the 
inverse proper radius 1/Ru [or l/{w{'Krc)Ru) for an arbitrary normalization of the warp 
factor.] 

Integrating out the universal extra dimensions, the Dirac neutrino mass matrix induced 
after electroweak symmetry breaking is 

= ^^ y^ (vrrj , (4.4) 

where Vh = 174 GeV is the Higgs VEV, and the denominator represents the square-root 
of the volume of the {T"^ /Z^ x [S^ / Z2) orbifold. As we will see in the next section, the 
neutrino mass eigenvalues are largely determined by the hierarchy among the ^^(Trrc), 
while the mixing angles are determined by the flavor structure of the couplings and 
hij. 



5 Zero-mode Profiles of Gauge-Singlet Fermions 

In this section we derive the profiles of the neutrino zero modes, which determine the 
neutrino mass matrix according to Eq. ([4.4|) . These depend on the VEV of and therefore 
our first task is to determine the solution to the equation of motion that follows from 
the A/'-independent part of Eqs. ( p.3|) and (|2.5| ). 

5.1 The bulk VEV 

We will be interested in a region of parameter space where the (p field VEV varies slowly 
in the bulk of the 7th dimension (with the exception of two narrow regions close to 
the branes), so that to a good approximation it simply gives a contribution to the bulk 
cosmological constant. Thus, we use the explicit form for the warp factor, Eq. ( |2.2| ). 

The negative mass-squared of if implies that a nonzero VEV for is energetically 
favored, but at the same time is an odd field under the orbifold identification, and 
therefore its VEV must satisfy the boundary conditions 

(y.(0)) = (^(vrrj) = . (5.1) 
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In terms of the rescaled VEV, 



which has mass dimension +1, the equation of motion is 

0_6.|^A,/3-M^/. (5.3) 

subject to the boundary conditions /(O) = /(vrrc) = 0. We assume A; > 0. In Eq. 
we have neglected possible higher dimension operators. This is justified as long as the 
effective field theory description is valid, that is as long as k and are well below M^,. 

Eq. ( p.3| ) describes the motion of a particle in the potential 

V{f) = -^f + -fP , (5.4) 

in the presence of an antz-friction term proportional to k. Thus, we are looking for 
trajectories in which the particle starts at the bottom of the potential / = with some 
initial velocity, climbs the potential up to a certain point and then rolls down back to 
/ = 0. The anti- friction term puts energy into the system, so it is conceivable that for a 
sufficiently large /c, no matter how small the initial velocity, the particle will gain enough 
energy to overcome the potential barrier; in this case the only solution that satisfies the 
boundary conditions is the trivial one f{z) = 0. 

We first determine the restrictions in parameter space for nontrivial solutions to exist. 
In order to do this it is convenient to set f{z) = ^J\, e~'^^^ f{z), so that Eq. ( p.3|) 
becomes 

g = My^^~f - [Ml - 9P) / . (5.5) 

describing now frictionless motion in a potential whose slope decreases with time. For 
k > M^/3, the "motion" starts at 2; = from the maximum / = of a continuosly 
decreasing potential, so that f{z) = is the only solution satisfying /(vrrc) = 0. Therefore, 
Eq. (|5.5| ) can have a nontrivial solution which satisfies /(O) = /(vrrc) = only if 

k<^ . (5.6) 

This necessary condition is not sufficient for the existence of solutions with /(O) = 
f{TTrc) = 0. Another necessary condition can be derived as follows. If we neglect the 
first term on the right-hand-side of Eq. (p. 51) the fictitious particle feels just a harmonic 
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f 




Figure 1: Scalar profile when M^p/k = 10, irrcM^ = 70. The plateau is at / M^/y^^. The 
points Zl and Z2 define the boundaries between regions I, II and III. In the main text we consider 
the case where zi/z2 <^ 1. 

oscillator potential and returns to the origin after a "time" z = i^{M'^ — 9A;^)~^/^. The ef- 
fect of the neglected term is always to increase the oscillation period. Hence, the boundary 
condition /(vrrc) = requires 

1 < - 9fc^ (5.7) 

c 

In the Appendix we prove that the conditions ( |5.6| ) and ( |5.7| ) are also sufficient for the 
existence of nontrivial solutions. We conclude that the VEV {^p{z)) is non-zero in the 
interval < z < ttVc for a substantial range of values of and k. 

For the special case of a flat and large z-dimension, k ^ l/(7rrc) ^ M^, the solution 
is given approximately by |31| 



f{z) ^ —j== tanh tanh ^"^^^^ 
V 2 V 2 



(5,i 



This solution is essentially constant except in the region of size ~ around the end- 
points. We will in general be interested in the parameter range in which k is not negligible, 
but where l/(7rrc) <C M^. The latter hierarchy will need to be only one to two orders of 
magnitude to explain the smallness of the neutrino masses relative to the weak scale. As 
we show in the Appendix, the solutions in this case are qualitatively similar to the flat 
case k = 0. We show a typical numerical solution in Fig. ||. 
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5.2 Bulk fermions 



2 J j ^3 



- hA,f - -j^f" e° , (5.9) 



In the presence of the (p VEV, the A/" fields have non-trivial profiles along the z dimension. 
The three Af zero-modes defined in Eq. (|4.1|) are a solution to the set of equations = 
1, 2, 3, with i fixed and j summed over): 

where we have again neglected possible higher dimension operators in Eq. ( |2.5| ). We also 
set /" = dPf/dz^ where / was defined in Eq. (|5.2|). It will be useful to factor out the 
leading order solution in powers of M^/M^,, namely the solution in the absence of the last 
term in Eq. (|5.9|) , by defining new functions Ci{z) through 

where 

S{z)^ fdUiC) . (5.11) 
The Cj's satisfy the following differential equations 

dz ~ Uf ^ ' 
which we now solve in the limit tiTcM^ ^ 1. 

Given the general features of the / profile discussed in the previous subsection, it is 
convenient to separate the analysis in the three regions shown in Fig. |l|. We first note 
that in region II, /" is exponentially small and therefore all the Cj's remain essentially 
constant throughout it. The differential equations ( ^.12|) are nontrivial in regions I and 
III. In region I the integral Eq. (|5.11| ) is of order A^^/^. In region III this integral is much 
larger, so that the important features of the Ci{z) are determined in this region as follows: 
for Z2 < z < TCTc the integral expression for S{z), Eq. (|5.11|) , can be replaced, to a good 
approximation, by S{Tcrc), which is itself of order A~^/^7rrcM^. Therefore, in this region 
there will be an exponential hierarchy among the various terms on the right-hand-side of 
Eqs. ( |5.12| ), provided /ijA~^/^7rrcM^ ^ 1. Without loss of generality we can assume the 
ordering hi > h2 > > 0. If we keep only the leading terms, then Eqs. ( [5.12[ ) in region 
III reduce to 



dCs ^ ^33 J// 

dz ~ M2 
dci h 



f"c3 (5.13) 



dz M? 



i^fn^{K-H)s(.r,) fori = 1,2. (5.14) 
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We first solve for ^3. It will be sufficient to work to zero-th order in M^p/M^, so that 
from Eq. ( ^.13|) we have 03(2;) = const. Imposing the normalization condition Eq. 



and in the limit l/(7rrc) <^ k < where we can evaluate S{z) by setting f{z) 

M^pl J \p throughout the region of integration, we find from Eq. ( |5.10| ) that 



^"(Trre) ^ Wyrr, (2M^ - k\ e-^^^--'^'^) , (5.15) 



where only one combination of the parameters in the Lagrangian for Lp appears, 

= ^j=^ . (5.16) 



We will see in the next section that the expression Eq. ( p.l5| ) will lead to exponential 
suppression of neutrino masses provided only that M<^ > k/2. 

Naively, one might think that for the other two generations there will be a greater 
exponential suppression controlled by their larger Yukawa couplings hi 2- We now show 
that this is not the case, due to the presence of the higher dimension operators that couple 
A/"^'^ to A/"^ in Eq. ( |2.5| ). To see this we need to solve for ci,2{z) in region III to first order 
in M^/M^. By using the zero-th order solution for 03(2;) ~ c^IttVc) in Eq. ( |5.14| ) we obtain 



c,{z) ^ Ci{z2) - C3(7rr,)^e('^-'^3)5K.) j^/^^) _ ^/(^^^j ^ ^ ^ ^^2, (5.17) 

from which we find 

^sl^rrc) C3(7rrJ 

But the first term here is exponentially small compared to the second term. Furthermore, 
f'{z2) is exponentially small compared to /'{nrc) and we finally obtain 

^iil^r^-I^rUr) (5 19) 

Because f{z) varies from 0{\-'^l'^M^) to zero over a distance of order 1/M^, it follows 
that —/'(nrc) = 0{X~^^'^M^). This provides the promised result. The zero-mode wave 
functions ^1 2(^'"c) are suppressed compared to ^^(Trrc), but only by a quantity 



2 



^--^-O^] . (5.20) 
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We note that in Eq. (|5.18| ) the leading term in 1/M* is exponentially suppressed while 
the sub leading term is not, so that it dominates as long as his is not extremely small 
compared with unity. The reader might wonder whether this signals a breakdown in the 
effective theory description, which relies on the convergence of the expansion in 1/M^,. 
There is no reason to worry: the term suppressed by in Eq. ( |5.18| ) comes actually from 
the leading flavor off-diagonal operator. All other terms suppressed by higher powers of 
give just small corrections to e. 



6 Neutrino Oscillations 

We are now equipped with all the tools necessary for analyzing the neutrino mass spectrum 
and the ensuing neutrino oscillations. 



6.1 Neutrino mass matrix 

As we discussed in the introduction (Section 1), the neutrino masses are of the Dirac 
type, the Majorana masses being forbidden by the symmetry under rotations of the two 
universal extra dimensions. Below the electroweak scale, the effective four- dimensional 
theory contains three left-handed neutrinos, z/|^, which are the neutral zero-modes of C[^, 
and three right-handed neutrinos, N}^, which are the zero-modes of A/'^'^'* with respect to 
the two universal dimensions. In the weak eigenstate basis. 



ulM.Nr + h.c. 



(6.1) 



the neutrino mass matrix derived in Eqs. (|4.4| ), ( ^.15 ) and (|5.19|) is given by 
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mo 



\13 \ 



e "'23'^ AT 



\23 



(6.2) 



where e -C 1 is defined in Eq. 



mo = Vhe 



and the scale of the neutrino masses is set by 

1/2 



-Trrcl M^-fc/2) 



/ 2M^ - k 



(6.3) 



.TxRuMj \ M, J 

As we will see in Section 7, the factors in parenthesis are expected to provide a suppression 
of no more than one to two orders of magnitude, so that the neutrino mass scale is 
explained by the first exponential. 
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The neutrino mass matrix is diagonalized by unitary transformations: 



uIm^Ur = diag(m^,^, m^^, m^^] 



(6.4) 



The unitary matrix describing neutrino oscillations, Ul, and the physical neutrino masses, 
m^., may be found by expanding in powers of e = 0{\~^^'^M^/M^). The largest physical 
neutrino squared-mass is 



= ml 



A 



+ 



A 



+ 



A 



It is convenient to use the following identity, valid up to corrections of order e^: 



-(/23/O* ^13 \ 



V 1^13 P + 1^33 P ^23 



V -h -{hshy I33 J 
where we use the following notation 



( eniii emi2 \ / 1 

£77121 £'"^22 

V / 



1 -e02 
V e^i e^2 1 J 



/'l \13|2 I I \23|2 _|_ I \ 



33|2 



1/2 



i,j = 1,2,3 



\i3 
A 



(|A)^r + |A 



33|2 
'AT I 



1/2 



1,3, 



(6.5) 



(6.6) 



3 ' 



1,2 . 



(6.7) 



The first matrix that appears in Eq. ( |6.6D can be shown to be unitary, while the last 
matrix is unitary to leading order in e^. The elements of the block-diagonal mass matrix 
shown in Eq. ( |6.6|) are given by 



mii 



"^!/3^i3 (^li^3 — ^3j^l 



h (hshi — hshi) + ^3 (^33^2* — hshi] 



, 2 = 1,2. 



(6.^ 



Eq. ( |6.6|) shows that the m^^ and m^^ physical neutrino masses are of order emu^, 
and generically are non-degenerate. We choose m^-^ < m^^. These can be computed 
straightforwardly by diagonalizing m (the 2x2 matrix whose elements are given by rhij): 



V^diag 



(6.9) 
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where V is a unitary 2x2 matrix. Its Vij elements depend only on ratios of A)^'s and 
on \hi3/h23\- This dependence can be computed straightforwardly using Eq. (|6.8|), but is 
cumbersome and we do not display it here. 

The unitary 3x3 matrix Ul is then given by the product of the first matrix on the 
right-hand-side of Eq. ( |6.6|) with 

/ V f] \ 

(6.10) 




The third column entries of Ul, 



[1 + 0{e')] 



13|2 , |X23|2 , |3^33|2\V2 



1,2 



(6.11) 



are relevant for atmospheric neutrino oscillations, as discussed below. The first two 
columns of Ul have entries that also depend (to leading order in e) only on ratios of 
A^'s and on \hi3/h23\. The ratio of the U}} and entries, which is relevant for solar 
neutrino oscillations, is given by 

tt12 r* \r 1 1 

(6.12) 



This ratio is typically of order unity if most of the A^'s have the same order of magnitude. 
In particular, this is true for the phenomenologically favored case discussed below, which 
is near the bi- maximal mixing of three neutrinos [^ . 







^21^3 — ViiZi/23 


UF 




Villi, — V21/1/23 



6.2 Predictions and experimental constraints 

The neutrino oscillation data constrains the differences of neutrino squared-masses. With 
three neutrinos, the atmospheric and solar neutrino data can be fit nicely (the LSND data 
cannot be accomodated). 

The atmospheric neutrino data require a larger mass splitting compared to the solar 
data so that m^^ is determined to leading order in e^. The range for the mass-square 
difference (Am^)atm given by the global fit to the data obtained in atmospheric 
neutrino experiments and in the CHOOZ reactor experiment is 

ml^ ^ {Am^)Zl = 1.5 - 6.0 x lO^^ eV^ at 99% CL . (6.13) 

As long as the combination of Yukawa couplings shown in Eq. (|6.5| ) is not smaller than 
unity by many orders of magnitude, Eq. ( |6.3|) requires 



"■'-2 



30 



(6.14) 
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where we anticipate that the terms in parenthesis in Eq. (|6.3|) are of order unity. 

The solar neutrino oscillations are controlled by the mass-square difference of the 
lighter neutrinos: 

(Am^)soiar = - ml^ ■ (6.15) 

Based on the reasonable assumptions that |/i23|, |^23| = and that most Yukawa 

couplings X'^ have the same order of magnitude, we find the following prediction for the 
solar neutrino oscillation scale: 



(Am^) 



'atm 



0{e') , (6.16) 



where e is given by Eq. ( p.2U| ). In the absence of fine-tuning, the (p scalar mass-squared 
can be no smaller than the one-loop contribution in the seven- dimensional theory: SM^ ~ 
(9[A^M^/(1287r^)]. Also, we expect that the effective theory description starts breaking 
down when approaches M^,. This leads to an allowed range for e [see Eq. ( p.20| )]: 

Oho-'^JX^ < t < O(l0-^/J\) , (6.17) 



with the values near the upper end being preferred if the mass parameter is not too 
much smaller than the fundamental scale M*. If A,^ is of order unity, the ratio shown in 
Eq. ( |6.16| ) is then of order 10~^. The generic prediction of our model for the solar neutrino 
scale then becomes 

(Am2),„i,, ^ 10-^ - 10-^ eV^ . (6.18) 

This prediction fits well the range currently allowed for the large-mixing angle MSW 
solution to the solar neutrino problem |^4|. The large uncertainties due to the unknown 
values of various parameters do not allow us to rule out the long oscillation wavelength 
solution, whose fit to the data prefers (Am^)soiar ~ 10~^ eV^. On the other hand, the 
vacuum oscillation solution requires (Am^)soiar ~ 10^^° eV^, so that it is disfavored within 
our model. 

The elements of the unitary matrix U l are constrained by the solar, atmospheric and 
reactor neutrino experiments. The element U]^ is most tightly constrained by the global 
fit to the CHOOZ and solar neutrino data, which gives {Uff < 6.5 x 10"^ at 99% CL 
35| . Applied to our case, this translates into a mild restriction: 



\23|2 I I \33|2 



< 0.26 . (6.19) 
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The atmospheric neutrino data strongly favors pure z/^ ^ oscillations with a mixing 
angle satisfying sin^ 2^23 > 0.83 at the 99% confidence level [^]. The expression of this 
mixing angle in terms of the Yukawa couplings is the same in our model as in the five- 



dimensional model of Grossman and Neubert |T^, and gives: 



0.65 < 



\23 



\33 



< 1.55 



(6.20) 



This is not a particularly strong constraint, given that the Yukawa couplings are typically 
expected to be of the same order of magnitude. However, if the march of the atmospheric 
neutrino data towards maximal ^ Vj. mixing continues, then it will be necessary to 
find a more detailed explanation for why the ratio of Yukawa couplings shown above is so 
close to unity. 

Finally, the global fits of the large mixing angle MSW solution to the solar neutrino 
problem require 



0.2 < 



< 



(6.21) 



This constraint is naturally accomodated by order unity values of the parameters entering 
the first two columns of the matrix shown in Eq. (|6.12|) . 

In summary, we can explain the hierarchy between the electroweak scale and the scale 
relevant for the atmospheric neutrino data by the first exponential factor in Eq. (|6.3| ). 
We can further explain the small hierarchy between this scale and the solar neutrino 
scale associated with the large-mixing angle MSW solution based on the typical size 
of the higher-dimension operators controlling neutrino flavor mixing. Furthermore, the 
currently allowed ranges for the mixing angles are natural if the various couplings of our 
model are of order unity. 



7 Mass Parameters and Warping 

Having observed that the smallness of the neutrino mass scale relative to the electroweak 
scale is explained by Eq. (|6.3| ) with the exponent given by Eq. ( |6.14|) , we now comment on 



the other mass scales of the model. These are the scale M^, at which the seven-dimensional 
theory, including its six- dimensional component, breaks down, the parameter k related to 
the seven-dimensional cosmological constant, the size of the universal extra dimensions, 
and the size of the warped dimension. We also introduced the mass of the bulk 
scalar that localizes the neutrino fields away from the standard model brane, and the 
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Higgs mass parameter Mh- These scales give rise to other, derived scales such as the 
electroweak scale Vh- 

There are several important relations among the above parameters. We assume a 
normalization of the warp factor as in Eq. ( p.2|) . (For a different normalization it is only 
necessary to interpret Ru in the following formulas as the proper radius of the universal 
extra dimensions at the standard model brane.) The first relation follows from the fact 
that the observed standard model gauge couplings, collectively denoted by g^, (as well as 
the top Yukawa coupling) are of order one. Writing a typical six-dimensional standard 
model gauge coupling in C^m of Eq. ( |2.3| ) as qq/M^, where is dimensionless, the observed 
four-dimensional gauge coupling is given by 

1 ~ f74 = (nR^M.y'ge . (7.22) 

If the standard model gauge interactions become strong at the scale M*, then ~ 1287r^ 



37|. In this case, the product RuM^, is of order vT287r. For a range of values of Qq > 1, 



will be somewhat above 1/Ru, providing a finite range of validity for the six- dimensional 
standard model. In the following we assume that C(l) < RuM^ ^ C(10).Q 

A second relation follows from naturalness considerations with respect to the Higgs 
mass. The Higgs mass parameter must be below M^, for the effective theory descrip- 
tion to be valid. However, on naturality grounds, Mh cannot be much smaller than the 
one loop corrections. They can be roughly estimated by cutting off the quartically diver- 
gent one-loop integrals, such as the one arising from the quartic Higgs self-interaction, at 
the breakdown scale of the effective, six- dimensional theory. This yields 




M,>Mh> 5Mh ~ \ TT^.M, , (7.23) 



where Ag is a dimensionless coupling in the six- dimensional theory. Using this estimate, 
it follows that the Higgs VEV is given by 



Vh 



T^RuM^X^^'^ 



Mh > -j^{RuM.fRz' . (7.24) 



^°An interpretation based on the AdS/CFT correspondence |Q may be useful if RuM^ < 1, in which 
case the present six-dimensional description is not applicable. We note that the Zs symmetry that lies at 
the heart of the remarkable proton decay suppression pointed out in and which also forbids Majorana 
masses, is expected to remain valid. Thus, even though it would not be possible to talk about a six- 
dimensional standard model, the six-dimensional structure still has observable effects in the low-energy 
four dimensional theory. 



24 



If RuM^ ^ then Eq. (|7.24|) can be read as an upper bound on in the TeV range. 
But the electroweak precision measurements impose a lower bound 1/ Ru ^ 0.5 TeV [|l|. 
It then follows from RuM^, ^ (9(10) that is in the TeV range. 

Finally we note that the parameters and k are expected to satisfy relations anal- 
ogous to Eq. (|7.23|) (though, as noted in subsection 5.1, we must have k < M^l?i). Our 
discussion of the solar and atmospheric neutrino data in subsection 6.1 indicated that 
must be approximately an order of magnitude below M^, which is consistent with a 
naturality estimate analogous to Eq. ( [7.23|) , if the corresponding coupling is of order 
unity or smaller. 

There are now various possibilities depending on how large the warping is. First 
suppose that the warping is no more than mild: e^"^"^ ~ 1. In this case, the weakness of 
gravity must be attributed to some suppression that lies beyond the seven-dimensional 
theory presented here, such as a few other flat dimensions accessible only to gravity along 
the lines of [Tl| . It is important to note that the mechanism for suppressing the neutrino 
masses presented here is independent of any such extension. From Eq. ( |6.14| ), we now 
have vrrcM^ 30. With roughly an order of magnitude below M^,, the inverse size of 
the seventh dimension, l/vrrc, is of order 10 GeV. 

A more interesting possibility is that -nkr^. ^ 1. In this case it is possible to explain the 



weakness of gravity a la Randall and Sundrum |jT2[ within the seven- dimensional model. 
There are, however, some differences arising from the existence of the universal extra 
dimensions. The four-dimensional Planck mass is now related to by 



Ak 



Ml, « ^"*^7^"^ e^---" . (7.25) 



Suppose that k is on the order of (but somewhat less than) M* . Then if 1 and are 
in the TeV range as expected from the previous considerations,!^ it follows from Eq. ( [7.25[ ) 
that 6^'^'"'=^ ~ 10^^, which translates into -nkrc ~ 20. (Note that there is an extra factor 
of two in the exponent compared to the five-dimensional Randall-Sundrum model.) We 
then see from Eq. ( |6.14|) that TrM^r^ ~ 40, and the inverse size of the seventh dimension, 
1 /Trrc, is again of order 10 GeV. The lightest spin-2 KK mode with momentum along the 

^^Note that we have chosen to measure all mass scales at the standard model brane with respect to 
the corresponding induced metric. Had we measured them with respect to a metric rescaled by e^^"^" as 
in Ref. [|l2j , we would have concluded that and the other "fundamental" parameters of the seven 
dimensional theory are of order Mp. The cutoff on the effective six dimensional standard model would, 
however, remain in the TeV range, being given by M,e~'^*"'''. None of the physical conclusions described 
here would change. 



25 



warped dimension has a mass of about 4/c, which is roughly of the same order as the mass 
of the first KK modes of the standard model fields. It will be interesting to see which KK 
modes will be discovered first if this model is realized in nature. 



8 Conclusions 

We have presented a higher- dimensional mechanism for generating a realistic neutrino 
mass spectrum. The smallness of the neutrino masses compared with the electroweak scale 
is explained by an exponential suppression of the right-handed neutrino wave functions 
on the standard-model brane. The hierarchy between the mass scales associated with the 
Ut- and u^, z/^ transitions, measured by the atmospheric and solar neutrino data, 
respectively, is limited by the effect of fiavor-non-diagonal, higher-dimension operators. 
As a result, the mass scale of the solar neutrino oscillations fits well the large-mixing- 
angle MSW solution. Furthermore, the neutrino mixing angles are naturally large if no 
large hierarchies between the neutrino Yukawa couplings occur. This is an important 
result in view of the fact that a majority of the models in the literature (for a recent 



review see |39|) can accomodate only the small- mixing- angle MSW or vacuum solutions 
to the solar neutrino problem, which are less favored by the data. In addition, the 
seemingly "maximal" mixing required by the atmosperic neutrino data is consistent with 
our mechanism for a reasonably large range of parameters. 

While our mechanism has been developed in the framework of the six-dimensional 
standard model, it is worth pointing out that it relies fundamentally on four ingredients 
that could naturally be present in a more general class of higher-dimensional theories: 
(i) Three right-handed neutrinos; 

{ii) A symmetry structure forbidding Majorana neutrino masses; 

(iii) A spatial dimension compactified on /Z2 and accessible to the right-handed neu- 
trinos but not to the standard model fields; 

[iv) A bulk (effective) scalar field which is odd under the Z2 orbifold transformation, has 
a VEV, and couples to the right-handed neutrinos. 

Ingredients {€) and {ii) are automatically present in the six-dimensional standard 
model, as required by the six-dimensional gravitational anomaly cancellation ^ and 
the Zg rotational symmetry of the two universal extra dimensions . We were then led to 
consider the six-dimensional standard model localized in a seventh dimension satisfying 

{iii) and {iv). 
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More generally, ingredient (ii) could be enforced by lepton number conservation in a 
variety of models, and its experimental test is the absence of neutrino-less double-beta 
decay. Ingredients {Hi) and (iv) could be present in 4+1 dimensional models, as suggested 
in Ref. ||31[ as a possible source of fermion mass hierarchies. We have generalized this 
construction by allowing a warped metric. We have analyzed the conjectured restriction 



on the number of gauge-singlet fermions |]T3|, and have found that in 6+1 dimensions 
the local and global gravitational anomalies cancel independently of the number of right- 
handed neutrinos, provided the fermion content is free of six-dimensional anomalies. (In 
4+1 dimensions no restrictions arise because there are no local or global gravitational 
anomalies in four or five dimensions.) 

Finally, it is worth noting that the model presented here has another intriguing feature. 
In the most appealing version of the model, the extra dimension that leads to an expo- 
nential suppression of the right-handed neutrino wave functions also solves the hierarchy 



problem along the lines of Randall and Sundrum ]T2[, while explaining proton stability 
based on six-dimensional Lorentz invariance as in Ref. 0. In this case collider searches 
at the TeV scale will reveal graviton Kaluza-Klein modes in addition to the Kaluza-Klein 
modes of the standard model fields. 
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Appendix 

In this Appendix we prove that the real scalar (p has a nonzero VEV whenever the 
parameters k, and satisfy the conditions ( ^.6|) and ( p.7|) . The rescaled VEV u{z) = 
^Jx^ f{z) must satisfy the z-independent equation 

This equation describes the mechanical motion in the potential —M^lu"^ — 1)^/4, and in 
the presence of an an^^-friction term proportional to k. The z coordinate plays the role 
of "time" variable. 
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We will prove the existence of a solution to ( |A.1|) , satisfying the boundary conditions 
m(0) = u{TTrc) = 0, by analyzing the "flows" in the equivalent first order system: 

du _ 
dz 

- QkP^ + Ml(u'-u) . (A.2) 



dz 

In the (m, Pu) plane, the solution sought corresponds to a flow that starts somewhere on 
the u = line and, after time vrrc, comes back to this line (because of the symmetry 
u —u, we can restrict u > 0.) More precisely, we shall prove that for any k < 
there are flows starting and ending on the u = line, with total elapsed times (vrrc) 
ranging from it /u to infinity, where 



u = J Ml - 9P . (A.3) 



The system ( |A.2| ) has two fixed points, (m, P„) = (0, 0) and (1, 0), near which it may be 



linearized. In the vicinity of (0, 0), Eqs. ( |A.2| ) has the following explicit solution satisfying 
the initial conditions {u,Pu) = (0, Pq > 0): 

uiz) = — e sma;^ , 

UJ 

P^[z) = u{z) (Sk + iutan~^ iuz'j . (A.4) 

Since u vanishes at z = n/u, the existence of solutions is established for a separation 
between the branes of ttTc h juj [see Eq. (|5.7|)], at least when Pq ^ such that the 
linear approximation is reliable. 

A similar linear analysis around (1,0) shows that there is an attractive fiow that 
approaches the fixed point (1,0) from the region m < 1, P„ > 0, and a repulsive fiow in 
the region m < 1, P„ < 0. For either one of them, the total time is infinite. This shows 
that adjacent fiows that come arbitrarily close to (1, 0) will spend an arbitrarily long time 
in its vicinity. Note that whenever the fiow enters the region u < 1, P„ < 0, it will be 
driven to the m = line, since the only possibility is to roll down the slope of the potential. 

We now prove the possibility of reaching these fiows from a point (0, P„ > 0), based on 
the following remarks. First, Eq. ( |A.2| ) shows that along the line P„ = 0, and for m < 1, 



the trajectories fiow vertically downward (see Fig.^. Second, the trajectories passing 
through the u = \ line fiow into the region u > \ whenever P„ > 0, and then they are 
driven to m ^ oo due to the negative slope of the potential. Tracing these trajectories 



28 



0.5- 



6k 



-0 . 5 




. 2 



. 4 



. 6 



. 8 



U 



Figure 2: The flow diagram when Mt^/'ik = 4. 



back in time, it follows from the first remark that they necessarilly cross the m = 0, 
-Pn > half-line. However, they cannot cross the origin because the flows cannot stop 
except at fixed points, and the point (0,0) cannot be reached [for small enough P„ > 
the trajectories stay only in the vicinity of (0,0), see Eqs. (|A.4|) ]. Note that the flows 
cannot cross due to the uniqueness of the solutions to ordinary differential equations such 
as Eq. (X]2). Third, tracing back in time the trajectory attracted to the fixed point (1, 0), 
the same argument as above shows that it crosses the line m = 0, P„ > 0. This "critical" 
trajectory, which starts at some critical (0, Pc > 0) and is attracted to (1,0), sets a 
boundary between qualitatively different flows. Finally, from any point (m,0), m < 1, the 
trajectory traced back in time crosses the line m = at some < P^ < Pc. Thus, a point 
(u, 0) with u arbitrarily close to 1 corresponds to a trajectory that spends an arbitrarily 
long time in the vicinity of (1, 0). 

We have shown so far that there are solutions for vrrc close to either ti / uj oi infinity. 
By continuity, it follows that there are solutions to Eq. (|A.1|) , satisfying the boundary 
conditions m(0) = uijiTc) = 0, for any > l/cj. This completes the proof that Eqs. ( p.6|) 
and ( |5.7|) are necessary and sufficient conditions for the existence of a nonzero (p VEV. 
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